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Low-mass, high-Q, silicon nitride nanostrings are at the cutting edge of nanomechanical devices for sensing 
applications. Here we show that the addition of a chemically functionalizable gold overlayer does not adversely 
affect the Q of the fundamental out-of-plane mode. Instead the device retains its mechanical responsiveness 
while gaining sensitivity to molecular bonding. Furthermore, differences in thermal expansion within the bilayer 
give rise to internal stresses that can be electrically controlled. In particular, an alternating current (AC) excites 
resonant motion of the nanostring. This AC fhermoelastic actuation is simple, robust, and provides an integrated 
approach to sensor actuation. 



Silicon nitride nanostrings^El are an exciting class of me- 
chanical resonator, exceptional for their ultra-high quality fac- 
tors (Qs)P their tolerance to irregularities in the fabrication 
process (since their resonant frequencies depend only on the 
device length) E' their harmonic normal mode spectrum, and 
the ease with which their vibrational motion can be detected 
optically (by virtue of the low absorption of silicon nitride 
at visible wavelengths) PEHI But in order for these devices to 
be used efficiently in room temperature molecular sensing 
applicationspE^fhey should have a surface that is chemically 
functionalizable. Silicon nitride, however, is largely inert. 

The most common workaround, routinely carried out with 
silicon and silicon nitride atomic force microscope (AFM) 
cantilevers, is to coat the resonator with gold. Then a 
molecule of choice^ can be affixed to t he g old overlayer, 
attached by way of a thiol intermedi ary! 1 4 I 1 5 I The addition 
of the gold has no adverse consequences for static, stress- 
response measurements,^ but it generally does for measure- 
ments that are dynamic in nature. For example in the can- 
tilever case, metallization causes significant dissipation in the 
resonator,^ lowering its Q and reducing its usefulness as a 
sensor — although this can be circumvented by depositing only 
in areas that do not become stressed upon actuation! 18 ! 19 ! 

Accordingly, our initial expectation was that application of 
a uniform gold layer to a silicon nitride nanostring would ruin 
the device's ultra-high Q, which is typically on the order of 
10 5 but can be as high as 7 x 10 6 (see Ref. 3 ). Instead, we 
show here that a 53 nm thick metallic layer on top of a 250 nm 
thick silicon nitride nanostring does not adversely affect the Q 
of the fundamental frequency, even though it nearly doubles 
the total mass of the string. (Our longest gold-covered device, 
at 210 /mi, has Q > 1.6 x 10 5 for the fundamental mode 
and a total mass of s» 840 pg.) It does, however, reduce the 
Qs associated with higher harmonics, so that the Q-versus- 
resonant-frequency behavior that was previously flat (up to 
even-odd effects across modes) now decays with frequency. 
This change is characteristic of a string system whose domi- 
nant dissipation mechanism is no longer localized at the an- 
chor points. 8 

Another important feature of the gold layer is that it ren- 
ders the string sensitive to temperature through the bima- 
terial effect (a differential stress between the gold and sil- 
icon nitride layers due to different coefficients of thermal 
expansion). While this could prove useful in tempera- 



ture sensingSSor provide an independent measurement of the 
resonator temperature, 24 our focus lies elsewhere. In partic- 
ular, we demonstrate that it is possible not only to vary the 
temperature of the device electrically via direct current (DC) 
ohmic heating but to actuate resonant motion via alternating 
current (AC) thermoelastic heating. Both the DC and AC 
effects have practical uses for sensors: the first allows for 
device regeneration (i.e., it provides a mechanism to desorb 
molecules, returning the device to a clean state), which facili- 
tates thermal studies of affixed molecules^HSaand, as recently 
shown, 24 can provide exceptional frequency stability in ultra- 
low-noise nanomechanical sensing; the second provides inte- 
grated actuation 30 32 and obviates the need for a piezoelectric 
buzzer or other external driving mechanism. 

Fabrication of the bilayer nanostrings is straightforward, 
since only one lithography step is necessary. We start with 
stoichiometric silicon nitride (250 nm) deposited onto sili- 
con dioxide (2 /im) on a silicon handle (Rogue Valley Mi- 
crodevices). This is subsequently sputtered with 10 nm of 
chromium, followed by 43 nm of gold. (We refer to the re- 
sulting device as bilayer nanostring, since the intermediate 
chromium layer is thin and serves only an adhesion function.) 
Standard optical lithography is then performed. This is fol- 
lowed by wet etching of the gold and then chromium, reac- 
tive ion etching of the silicon nitride, and device release with 
a buffered oxide etch of the silicon dioxide. The resulting 




FIG. 1. False color scanning electron microscope (SEM) image of 
the 210 fim long gold and silicon nitride (purple) bilayer nanostring 
(scale bar 20 /im). Inset is an overhead view (scale bar 2 /tm). The 
gold layer is a 2.18 /im wide strip atop the 2.75 /tm wide string. 
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strings are shown in Fig. [T] The gold is slightly over-etched, 
although this is not important since the string's normal modes 
are insensitive to the exact dimensions. This fault tolerance 
during fabrication is advantageous for bulk processing. 

Measurements are performed by optical interferometry us- 
ing a 632. 8 nm laser focused onto the nanostring, as described 
elsewhere,^ 3 ! inside an optical access vacuum chamber in 
which the temperature of the sample stage can be controlled 
and accurately determined. The optical power incident on 
the device is 13 /iW — below the onset of optical heating ef- 
fects. The interferometric signal is amplified 25-fold after de- 
tection and is analyzed using a high frequency lock-in am- 
plifier (Zurich Instruments HF2LI), which allows us to ex- 
tract the peak frequency and Q by fitting to the power spec- 
tral density function (PSD) (see supplementary material of 
Ref. |S). This analysis is appropriate regardless of the actu- 
ation method: thermomechanical, external piezoelectric, or 
integrated ohmic (thermoelastic). In addition, our optical sys- 
tem can be scanned with respect to the sample using closed- 
loop piezo stages, 33 providing spatial maps of the amplitude 
and phase output by the Zurich lock-in. Hence, we can unam- 
biguously distinguish between in-plane and out-of -plane mo- 
tion. 

We emphasize that while we have verified the following re- 
sults using numerous nanostrings of varying lengths, for con- 
sistency we present data for a single 210 /im long by 2.75 fim 
wide bilayer nanostring with an initial resistance of 176 fl. 
The Qs of this bilayer nanostring in vacuum (w 2 x 10~ 6 torr) 
are shown in Fig. [2] These Qs are compared with out-of-plane 
modes of uncoated strings of a very similar geometry (red 
squares) from Ref. [8] It is immediately clear that the Q of the 
fundamental out-of-plane mode is essentially unchanged (al- 
though the resonant frequency itself is reduced by about 25% 
because of the addition of the gold). This device is therefore 
an excellent candidate for a functionalizable, high-Q nanome- 
chanical sensor. 

On the other hand, the higher order modes are damped 
as compared to the bare string. Both bare and gold bilayer 
devices demonstrate the alternating even/odd behavior previ- 
ously seen for silicon nitride membranes 34 and nanostrings, 8 
which has been reported to originate from dissipation due to 
phonon tunneling. 34 Ignoring this oscillation, the Q of the bare 
device is essentially flat with respect to frequency (or, equiv- 
alently, mode number). The Q of the gold bilayer strings, 
however, decays with frequency. A key difference appears to 
be the high thermal conductivity of the metallic layer, which 
allows for efficient transport of energy from within the string 
(dissipated as a result of internal stresses during bulk bending) 
to the outside. 

As we have previously shown, 8 the Q varies by mode ac- 
cording to 
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The tilde decorated quantities are damping coefficients from 
which the leading order frequency dependence has been fac- 
tored out. The viscous damping term, ^ Vlsc , can be ne- 
glected since the experimental chamber has been evacuated 
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FIG. 2. Q as a function of mode number for the 210 ^m long gold 
bilayer nanostring (orange circles are out-of-plane and blue triangles 
are in-plane) and a bare string of similar geometry (red squares) from 
Ref. [8] Here the string is actuated using the external piezoelectric. 
The blue and orange curves are fits to Eq. {T} with 7 VISC = 0. 



far past the point where the measured Q exhibits a pres- 
sure dependence. 4 In the case of the bare string, ^ anchor is 
the overwhelmingly dominant contribution; 8 energy is lost 
primarily in dissipative processes in the vicinity of the an- 
chor points, and the Q has no frequency dependence other 
than the even/odd mode effect. The Q of the gold bilayer 
string, however, shows a decaying tail that can only be fit 
with a non-vanishing 7 bulk term, attributable to bulk dissipa- 
tion. Specifically, we find from the out-of-plane data in Fig. 2 
that Wl 7 bulk / 7 anchor = 0.38 ± 0.08. Determination of the mi- 
croscopic bulk dissipation mechanism would require further 
detailed studies. 35 We remark that the Qs of the higher order 
modes are still rather impressive (> 10 4 ), suggesting that this 
device can be used for multimode sensingPEZl 

Beyond altering the dissipation, the gold overlayer adds 
functionality to the nanostring system by introducing bima- 
terial temperature dependence. In particular, we have ob- 
served that the resonant frequencies of the nanostring shift 
with heating. There are two scenarios to consider: (i) the 
sample stage is held in thermal equilibrium at a uniform tem- 
perature T = T room + AT; (ii) heating is highly localized in 
the string, and the rest of the system — including the supports 
and substrate — can be viewed as a heat bath at T room . The sec- 
ond scenario corresponds to the Joule heating case in which a 
current is passed through the metallic overlayer. 

In the first scenario, the frequencies have the form 
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where a k $ — a k — &o are the relative thermal expansion co- 
efficients measured with respect to that of the substrate, ao, 
and the parameters E k , A k are the corresponding Young's 
moduli and cross-sectional areas. (See Ref.l43land references 
therein.) These material parameters are defined for each layer 
with the indices k = 1 (silicon nitride), k = 2 (chromium), 
and k — 3 (gold). 

For small temperature changes, it is sufficient to consider 
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FIG. 3. (a) Relative frequency shift of the fundamental mode in vac- 
uum (actuated with external piezoelectric) versus temperature as the 
entire sample stage is heated, (b) Fundamental out-of-plane mode ac- 
tuated by the AC ohmic thermal effect, (c) Absolute frequency shift 
of the fundamental mode in vacuum versus applied voltage using AC 
ohmic actuation (red squares) and corresponding device resistance 
(blue circles). Fit (dashed grey curve) is to Eq. |6}. 
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We observe that the prefactor to AT is either positive or neg- 
ative depending on the value of an. For our device, the fre- 
quency of the resonator shifts upward with heating since the 
thermal expansion coefficient of the substrate, ao = 2.59 x 
1(T 6 K-\ is larger than (Y^ =1 £ fc A fe a fe )/(ELi E M = 
1.84 x 10 _6 K _1 , an effective alpha of the string that de- 
pends on the expansion coefficients of the constituent ma- 
terials, averaged layer-wise according to the product of the 
cross-sectional area and Young's modulus. Hence, the string 
is stretched tighter as the temperatures goes up; the relative ex- 
pansion of the string increases the effective tensile stress, and 
the resonances shift to higher frequencies. Putting reasonable 
estimates of the material parameters 38 into Eq. (|3), we predict 
1.30 x 1CT 4 K _1 for the slope in Fig .(3^a). This is consistent 
with the observed value, 1.10 x 10~ 4 K _1 . 

In the second scenario, resistance to an electrical current 
through the metallic layer of the string causes local heating. 
The substrate is unaltered and only the string deforms. The 
frequency shift then has the same form as Eq. |2) but with 
each dk,o replaced by its bare value a k and, if the current 
is alternating, with AT understood to be the time-averaged 
temperature variation. 

We can model the resistive heating in a simple way. If a 
constant power P is applied to the string, the total excess heat 



energy U comes into equilibrium when U = P — U/t± = 0. 
Here, r t j, represents the half-life for thermal energy to leak out 
of the string. The energy outflow can be computed as a sum 
of thermal currents (of the form j = — kVT) driven by a gra- 
dient VT w AT/(T/2). We make use of the device's overall 
heat capacity C = T£Li A kPk c k = 3.70 x 10~ 10 J/K to 
relate the change in heat energy to a corresponding change in 
temperature and thus to produce an estimate of 
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(Values of the density p k , specific heat c k , and thermal con- 
ductivity Kk are tabulated per layer in the supplementary 
materialPO) Note that r t h scales as L? , since the equilibra- 
tion process is essentially diffusive; this implies that, near 
resonance, ut^ oc L. Because the layers themselves come 
into local thermal equilibrium on a fast time scale of around 
Ct\jLw\K\ ~ 5ns (where t\ and w\ are the thickness and 
width of the silicon nitride respectively), we proceed as if the 
silicon nitride, chromium, and gold layers are at the same tem- 
perature instantaneously at each point along the string. This 
is a safe assumption up to a frequency scale l/(27r • 5 ns) « 
32 MHz, which is well above the oscillation frequency of any 
of the harmonic modes that we can measure in our device. 

The simple analysis we have outlined is primarily limited 
by the fact that we have supposed the temperature profile 
AT(x, t) = AT(i) to have no meaningful variation along the 
length of the string. Still, it reliably captures the fact that the 
thermal time scale is on the order of microseconds (rather than 
nanoseconds, as for the device described in Ref.[30l. orders of 
magnitude faster than the time for the mechanical energy to 
be transferred to the environment, as dictated by the high Q of 
the string: Q/f x 200 ms. 

In our AC local heating scenario, the power pumped into the 
string is given by P = V 2 / R, where R is the string's electri- 
cal resistance and V(t) — Vq coswt is the alternating voltage 
applied across the string. Since the silicon nitride is insulating 
(i?i « oo), we can estimate the resistance from the geometry 
of the device using standard values f2, r% for the resistivity of 
CrandAu: R = (£Li l/Rk)- 1 w L(A 2 /r 2 + A 3 /r 3 )- 1 = 
49 fL The 176 ft we measure experimentally is higher, mainly 
because it includes the resistance of the metallic supports. Ac- 
cordingly, the heat content of the string (ignoring the initial 
transient) varies in time as 
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with a phase lag tj> — tan -1 (2wr t h) that depends only on 
the thermal outflow time r t h. We can decompose the cor- 
responding temperature variation into a time-invariant back- 
ground (AT) ave = V 2 T t h/2Ci? over which is superposed the 
weaker, alternating behaviour (AT) a i t . Hence, the analogue 
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FIG. 4. (a) Frequency response of the first out-of-plane (lower 
frequency) and in-plane (higher frequency) modes in air using AC 
thermoelastic actuation (400 mV drive) with corresponding PSD fits 
(black) to extract the Q. (b) SEM of the gold bilayer nanostring after 
current induced failure at w 450 mV AC drive. Enough gold has 
come off the device to break the current path, resulting in insulating 
behavior. Scale bar is 2 /xm. 
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The relative frequency shift is quadratic in Vq and always 
negative — consistent with our observations in Fig. |3jc). In 
contrast to the relative frequency shift described by Eq. 
which is insensitive to the device length, the effect described 
in Eq. |6]) is predicted to scale as t±/C ~ L (and thus can be 
made arbitrary large by choosing a long enough string). The 
predicted coefficient, 6.82 V~ 2 , is about ten times larger than 
the value 0.648(9) V~ 2 obtained from a fit to the data. Theory 
and experiment can be made to coincide if instead we take the 
thermal time to be t± = 72 /is. The discrepancy likely stems 
from the over-simplified, uniform temperature profile along 
the length of the string. 

In the vicinity of the fundamental mode of our device, 
2wiT t h = 47r/iT t h = 768 3> 1, so that <fi has essentially 
achieved its maximum value, 7r/2. In this limit, (AT) a i, = 
(Vq/2CRuj) sin(2tdf), and the magnitude of the alternating 
component of the temperature with respect to the average goes 
as 1/wTth <C 1. What's astonishing is that this tiny alternat- 
ing local heating is sufficient to actuate the device and to do so 
very efficiently. Specifically, we are able to use this local heat- 
ing to actuate the bilayer nanostring by passing an AC current 
along the length of the device. Note that this is subtly differ- 
ent from the bimaterial effect in cantilevers and non-uniformly 
coated beams, where the gold exerts a force on the cantilever 
that pushes it towards the substrate.EEH i n sucn a scenario, 
the force is directional, and hence the AC drive voltage must 
be applied at frequency f/2 in order to produce heating and 
bending at /. In contrast, in our system there is no preference 
towards or away from the substrate as the gold expands, since 
it forms a continuous layer over the string. Heating and bend- 
ing occur twice per string oscillation, and therefore the applied 



drive voltage and resulting motion have the same frequency. 
We are also able to actuate the device at f/2, although w 1000 
times less efficiently — consistent with the above scenario. 

In Fig. [3jb), we show frequency sweeps using the refer- 
ence output of the Zurich lock-in amplifier to resonantly drive 
the nanostring. This integrated electronic actuation scheme 
is simple, compact, and robust. With root-mean-square AC 
voltages as small as 100 fj,V we can clearly detect resonant 
motion, while the maximum signal-to-noise is achieved with 
approximately 15 to 20 mV. At higher drives the nanostring 
becomes nonlinear. Yet even at quite high drives, in vacuum, 
we see no change in the resistance of the gold layer. Hence 
the local heating of the bilayer nanostring is nondestructive 
to the gold layer in the linear actuation regime and beyond. 
In addition, we note that the resonance frequency is entirely 
non-hysteretic at these actuation voltages — further evidence 
that the AC ohmic actuation scheme is robust. Only at volt- 
ages > 180 mV, approximately 10 times that required for op- 
timal signal-to-noise, is there an onset of resistance changes 
in the gold layer. The device continues to function with ap- 
plied voltages in excess of 350 mV, after which it rapidly be- 
comes an open circuit. Using the extracted value T,h = 72 [is, 
we can estimate that the upturn in resistance, appearing at 
Vq rj 180 mV in Fig. 3(c), begins when the device has been 
heated an additional AT sa 56 K and that burn out occurs 
when V « 450 mV and AT w 350 K. SEM of our device, 
shown in Fig. |4jb), reveals that the gold has only been re- 
moved from a small patch on the string surface yet has inter- 
rupted the current path. 

Finally we note that the AC ohmic actuation technique 
works for both higher-order harmonics in vacuum (not shown) 
as well as for the fundamental modes in air. In the ambi- 
ent pressure case, the nanostrings experience strong viscous 
damping, especially when the motion is perpendicular to the 
string's substantial width. The Q of the first out-of-plane 
mode is only 40 while that of the in-plane mode is 160 [see 
Fig.Qa)]. There is a greater cross-section for molecular colli- 
sions, and thus energy loss, in the out-of-plane direction. The 
voltages required to actuate the nanostring in air (« 400 mV) 
result in resistance changes and eventual failure. Improving 
the Q in viscous environment*^ 4 ^ would decrease the re- 
quired actuation voltages and prevent such failure. 

In conclusion, gold bilayer nanostrings are an enhancement 
to the already exciting system of silicon nitride nanostrings. 
They enable chemical functionalization and molecular detec- 
tion, provide a scheme for local heating to study thermal re- 
sponse of molecules or desorb unwanted molecules, and pro- 
vide simple and integrated actuation. And unlike traditional 
micromechanical resonatorsP^l they maintain their high qual- 
ity factors even with a continuous metallic coating. 
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SUPPLEMENTARY MATERIAL 

Device geometry 

The device is shown in cross section in Fig. |S-l| Each of the 
silicon nitride (k = 1), chromium (k = 2), and gold (k = 3) 
layers is of constant thickness tk and has an area Ak and a 
bending moment Ik given by 



A k = / dCw(0, 
J(k-i 

JCk-i 



(S.l) 



The limits of integration are the cumulative layer thicknesses 
= Ya=i^i- Tri e bending axis of the composite string 
is fixed at a height ( — 136 nm = 0.54ii (measured from 
the bottom surface) that minimizes the total bending constant 



Hence, the Ik values themselves de- 



pend indirectly on the Young's modulus Ek and Poisson ratio 
Vk of the material in each layer. 

Based on our best understanding of the wet- and dry-etch 
steps in the fabrication process, we propose that the layered 



6 




4x£ 



cross section 



w 3 = 2.18|im 



Cr 



Si 3 N 4 


=250nm 

1 





u = lOnm 



W! = 2.75 [im 



FIG. S-l. Bilayer nanostring viewed in cross section. Scale magni- 
fied four-fold in the vertical direction. 



string has a width profile of the form 

f ^3 ~ 2y /tf^(C3 - Z? if C2 < ^ < C3 

w(z) = lw 3 2^1 - (C 2 - z) 2 ifCi<z<C 2 (S.2) 
[ wx if < z < Ci 

The values wi = 2.75 /im and W3 = 2.18 /im are estimated 
from a (top-down) scanning electron microscopy image of the 
device. We use a pixel-contrast histogram analysis to perform 
edge detection. Finally, the areas and bending moments have 
the following values: 



(A U A 2 , A 3 ) = (68.8,2.16,9.08) x l(T i4 m 2 , 
{h,h,h) = (36.6,3.08,19.4) x 10" 28 m 4 . 



(S3) 



The device geometry and dimensions we have indicated 
correspond to a device of total mass 841 pg and heat capac- 
ity 3.70 x 10- 10 J/K. 



Normal Modes 

The layered beam is described by a Lagrangian 

C = £ dx (J^(d t z) 2 - U[z(x,t)]\ (S.4) 

where M/L — Y2k=i P^Ak is the (uniform) mass per unit 
length and 

1 1 3 

U = -oA x {d x zf + - E k A k (d x z) 4 



k=l 



k=l k 



is the elastic energy density (including tension, elongation, 
and bending terms) associated with its transverse deforma- 
tion. U penalizes any deviation from the equilibrium posi- 
tion z(x,t) = 0. At this stage, we ignore thermal expan- 
sion effects and restrict ourselves to a model that satisfies 
U[z(x,t)} = U[z{L~x,t)} and U[z(x,t)} = U[-z(x,t)]. 
In other words, we assume that the system is invariant under 
the transformations z(x) — > z(L ~ x) and z(x) —> —z(x), 
even though the latter symmetry can in principle be broken by 
the layered structure. 

We follow the usual prescription 



5C 



SC 



dt 5[dtz{x,t)] 8z(x,t) 
and make use of the variational identities 

S 



(S.6) 



5z(x, t) 
S 

Sz(x, t) 
5 



dx (d x z) 2 
dx (d x z) 
dx{d 2 x zf 



Sz(x, t) 

to obtain the equation of motion 

3 



-2d 2 x z 
-\2{d x z) 2 d 2 x z 
2{d x zf 



(S.7) 



Md z t z = Td z x z 



S(d x z) 2 diz - Bd% 



(S.8) 



Here, r = LA\a is the tension in the composite string, and 
S = i^Li E kA k and B = LY, 3 k=1 E k I k /(l - v 2 ) are 
its stretching and bending constants. It is sufficient to treat 
(d x z) 2 d 2 z, the term that is nonlinear in z, at the mean field 
level. We do so by replacing (d x z) 2 by its spatially averaged 
value 



dx (d x i 



(S.9) 



Equation ( |S.8[ ) differs from a pure wave equation by terms 
that are proportional to the Young's moduli. The relevant sit- 
uation here is when these are sufficiently small that the be- 
havior of the system is almost perfectly string-like. When 
the system is driven in the nth pure string mode z(x,t) — 
a n sin(n7ra;/L)e IIi; " t , it oscillates at an angular frequency 



1 

M 



L 2 



-Sfi 



L 4 



-B 



(S.10) 



For weak nonlinearity, the mean square displacement (a 2 ) 
H is small compared to L 2 t/S. In that limit, 



■I n 



1 

4M 



L 2 



L 4 



-B 



(S.H) 



Accordingly, high order modes exhibit a small n 3 correction 
away from a perfectly harmonic spectrum: 



In 



n 

2L 



M AL 3 JVm' 



(S.12) 
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TABLE S-I. Material properties of the layered nanostring and substrate. Unless otherwise specified, reference values for the bulk material are 
used (drawn from CRC Handbook of Chemistry and Physics, 92nd Edition). 











thermal 


specific 




thermal 


layer (#) 


density 


Young's modulus Poisson ratio 


expansion 


heat 


resistivity conductivity 




P 


E 




a 


c 


r 


K 




(g/cm 3 ) 


(GPa) 




(xlO 6 K _1 


) (J/kg-K) 


(nO-m) 


(W/m-K) 


Au(3) 


19.30 


79 


0.44 


14.2 


129.1 


22.14 


318 


Cr (2) 


7.19 


279 


0.21 


4.9 


448.2 


125 


93.9 


S13N4 (1) 


3.043 


265 a 


O.2C0 


1.27 a 


700 


00 


30 


Si (0) 








2.593 









a determined by multimode calibration measurements 
b from Ref.|4T1 
c from Ref.1421 



Thermal Expansion 

Expansion of the layers due to heating leads to a correction 
that enters the elongation term of the elastic energy density: 

1 1 3 

U = -A^idxz) 2 + E k A k [{d x zf - 2a k AT] 2 



k=l 



lilY^M*) 2 - (S.13) 



k=l 



Note that we are not attempting to model interfacial shearing 
and peeling stresses. We just assume that the layers are held 
together by a sufficiently large through-thickness spring con- 
stant. Expanding Eq. ( |S.13| l to give 



U = const 



1 

2L 



(r-L^B t 4 Qt AT) {d x zf 



fe=l 



B{d 2 x z) 



, (S.14) 



we see that the effect of heating is simply to shift the tension. 
Applying this observation to Eq. ( |S.l l| l gives 



/„(AT)-/„(0) 

/n(0) 



1 



8ML(f n /n)- 



E k A k a k AT 



k=l 



la A x 



1 

— E k A k a k AT. 

(S.15) 

The last line above is correct to the extent that f n /n « f\. 

Material properties 

The properties of silicon nitride thin films are quite differ- 
ent from those of the bulk material. For example, the density 



and coefficient of thermal expansion are generally lower. So 
we have taken care to find values that are appropriate to our 
particular devices. 

For the purposes of calibration, we prepare a bare silicon ni- 
tride nanostring with no metallic coating. The dimensions of 
this device are 215 /im (length), 2.1 /im (width), and 250 nm 
(thickness). The frequencies (in MHz) of the first six reso- 
nant modes are determined to be 1.19508, 2.38389, 3.58222, 
4.77447, 5.98263, and 7.19362. In this case, Eq. ((STTT} re- 
duces to 



(2Lf n ) 2 



a 






hn 4 ?r 2 


P 




a 


aE 




hn - 


P 


P 



p 12(1 - v 2 ) 



(S.16) 



A best fit of this form to the data produces 



P 
E 

P 



= 2.622(3) x 10 b m 2 ls\ 



8.7(8) x 10 7 m 2 /s 2 , 



(S.17) 



and hence E/a = 330(30). We have no direct measurement 
that closes this set of relations, but if we assume the manu- 
facturer's quoted value of 0.8 GPa for the tensile stress, then 
E = 265 GPa and p = 3.05 g/cm 3 . These values are consis- 
tent with other measurements in the literatures^ 

Furthermore, we obtain the relative frequency shift for each 
of modes 1, 3, and 5. The measurement is performed over 
a dense grid of temperature values from 293 K to 320 K and 
produces a linear slope 



E(a Q - ai) 
2o 



= 0.000219(1) K" 



(S.18) 



that is consistent across all three modes. Taking a = 2.59 x 
10 -6 K _1 as the coefficient of thermal expansion for the sili- 
con substrate, 42 we find a\ — 1.27 x 10 -6 K _1 . This is very 
close to the value 1.23 x 10~ 6 K^ 1 reported in Ref.|43l 



